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The analytical continuation average spectrum method (ASM) and maximum entropy (MaxEnt) 
method are applied to the dynamic response of a noninteracting resonant level model within the 
framework of the Kubo formula for electric conductivity. The frequency dependent conductivity is 
inferred from the imaginary time current-current correlation function for a wide range of tempera- 
tures, gate voltages and spectral densities representing the leads, and compared with exact results. 
We find that the MaxEnt provides more accurate results compared to the ASM over the full spectral 
range. 



I. INTRODUCTION 

The computation of real time correlation functions 
in many-body quantum systems is challenging due to 
the exponential complexity of evaluating exact quantum 
dynamicsP^ This is exemplified by the well-known dy- 
namical sign problem common to real-time Monte-Carlo 
techniques. The sign problem can be avoided in imagi- 
nary time and Wick rotation may be used to recover all 
the real time information and excitation spectra. How- 
ever, since the imaginary time correlation function must 
be determined numerically via Quantum Monte Carlo 
(QMC)Pthe rotation becomes numerically unstable and 
is highly sensitive to statistical errorsP 

Several approaches have been developed in order to 
circumvent this problem, such as the maximum entropy 
(MaxEnt) method where the optimal fitting of a data 
is defined in a Baycsian manner, which in general de- 
scribes the competition between the x 2 goodness of the 
fit and an entropic prior £P3 The MaxEnt approach 
has been applied successfully to a number of physically 
interesting problems pHUl A second approach is based 
on the notion of averaging over a sequences of possi- 
ble solutions. Such approach is called stochastic analyt- 
ical continuation method or average spectrum methods 
(ASM) J2M3 Recent application of the ASM and MaxEnt 
approaches argued that the ASM should be superior to 
the Ma xEnt a t least in its ability to resolve sharp spectral 
features! 28 * 29 ^ Examples included the calculation of the 
dynamical density fluctuations in liquid para-hydrogen 
and orf/io-deuteriuirPSl as well as spin dynamics in anti 
ferromagnetic Heisenberg spin champs 

In this work we apply both methods to study the dy- 
namic response of the well known resonant level model 30 
and exhibit the dc conductivity extracted from these ap- 
proaches as well its frequency dependence. Our approach 
to describe the conductivity is quite different from that 
discussed recently in the literature) 3 -^ in that it is not 
limited to a specific choice of the Hamiltonian; it does 
however, apply to equilibrium situations only. To assess 



the accuracy of both analytic continuation methods, we 
compare the results to the exact solution, which predicts 
a smooth broad spectrum. We find the MaxEnt approach 
provides an overall good agreement with the exact solu- 
tion for the entire range of frequencies and is particu- 
larly accurate at low frequencies near the dc regime. On 
the other hand, the ASM yields sharp, narrowed spec- 
tral features as well as spurious concave domains which 
are noticeably different from the exact results. Moreover, 
its low frequency predictions are somewhat less accurate 
compared to those of the MaxEnt method. 

The paper is organized as follows: In Section |TT] we 
describe an analytic continuation approach for electrical 
conductivity and briefly review two techniques to perform 
the Wick rotation: The ASM and the MaxEnt method. 
In Section IIIII we describe our model Hamiltonian and 
the exact solution for the electrical conductivity within 
Kubo's linear response theory. Section [TV| is devoted to 
present the results obtained from the ASM and MaxEnt 
approaches. Conclusions are given in Section M 



II. ANALYTICAL CONTINUATION 

The objective is to obtain the electrical conductivity 
from Kubo's linear response theoryP^l 



<t k (uj) = ft dte 
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where ji — -rKp is the inverse temperature, / is the cur- 



rent operator defined below (see Eq. 14 1, and the la- 
bel (••■)« represents the Kubo transform of the current- 
current correlation function defined byP^l 
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(I(t)I(0)) K = — y d\(I(t - »A)J(0)>. (2) 

In the above, h is the reduced Planck's constant and (• • • ) 
denotes an ensemble average. 
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To obtain the electrical conductivity, one requires the 
calculation of the current-current correlation function. 
C 0) = (l(i)J(O)}, in real time: 



C(t) 
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where _ff is the Hamiltonian of the problem and Q the 
canonical partition function. C (i) is related to the fre- 
quency response function D (w) by a simple Fourier rela- 
tion: 



C (t) = - 1 - / e- tut D (lu) du). 
2n 
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The frequency response function, D(lu), can be written 
in terms of the real part of the electrical conductivity, 
<j K (u>): 



5R{a K («)}, (5) 



where for uj — > 0. the dc conductivity is given by 
(%(^0)} = 0): 



a K (Q) = £l?(0) 



(6) 



To obtain a (lu), one requires the direct calculation of 
C (t) , which is an extremely difficult task for a typical 
model Hamiltonian used to describe transport in confined 
systems. However, in comparison the real time case, the 
calculate the corresponding imaginary-time function by 
QMC techniques is less tediousP^ In imaginary time, the 
corresponding autocorrelation function is also related to 
D (lu) . This relation is achieved by performing the Wick 
rotation, replacing t —> —it and using the detailed bal- 
ance relation D (—uj) = e~^ hu D (lu): 

oo 



where < r < h/3 and the imaginary time autocorrela- 
tion function can be calculated by modifying the Hciscn- 
berg equation of motion in Eq. (|3| : 

C(-iT) = ~Tr{e-P"en*Ie-i"f}. (8) 

Much of the difficulty in calculating C (t) is now shifted 
to that of inverting the above relation to obtain D(oj) 
from C(—ir). We refer to two approaches developed to 
invert this relation: The MaxEnt method and the ASM. 
The MaxEnt method selects the solution which maxi- 
mizes the posterior probability, or the probability of the 
solution D(ui) given a data set C(—ir). Using Bayes' 
theore m, one can show that the posterior probability is 
given b} 134 | 3S | 



where S is the information entropy and \ 2 is the standard 
mean squared deviation from the data. In the present 
study, we only consider the standard L-curve method to 
determine a™ In this context we regard a as a regular- 
ization parameter controlling the degree of smoothness of 
the solution, and the entropy as the regularizing function. 
Its value is selected by constructing a plot of log[— S] vs. 
logx 2 . This curve has a characteristic L-shape, and the 
corner of the L, or the point of maximum curvature, cor- 
responds to the value of a which is the best compromise 
between fitting the data and obtaining a smooth solu- 
tion. Our experience with analytic continuation has lead 
us to the conclusion that this approach is the most ro- 
bust. Once the regularization parameter is determined, 
we apply the approach described in Ref. |7] to obtain the 
statistically rigorous and unique fit for D(lu). 

The basic idea behind the ASM is to pick the final solu- 
tion for D(ui) as the average spectral function obtained by 
averaging over a posterior probability, V (D(lo)\C (—ir)) , 
instead of taking the value that maximizes this distribu- 
tion. Thus, fluctuations of the solution are allowed in the 
ASM: 



D(u) 
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The averaging is performed by a Monte Carlo procedure 
preserving known sum rules pSl Readers who are inter- 
ested in a more comprehensive discussion of the ASM 
are referred to Refs. l37l and l38l and to the specific imple- 
mentation of Ref. [Ml 



III. MODEL 

We consider the resonant Level model, which con- 
sists of a single quantum dot state coupled to two leads 
(fermionic baths). This model has been used extensively 
in understanding transport properties of non-interacting 
systems, most recently in a scmiclassical study of trans- 
port developed by Swenson et al. 3 - Since an exact solu- 
tion of the frequency dependent conductivity is available, 
this model is ideal for assessing the accuracy of the pro- 
posed analytic continuation approach. The Hamiltonian 
is given by (more details can be found in Ref. IBT)]) 



H 



d)d+ 2_ j e k c\c k + V t k (d^c k + h.c.) , (11) 
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V(D(oj)\C(-it)) cx exp(aS' - X 2 /2). 



(9) 



where d (d^) is the destruction (creation) operator of an 
electron on the dot, Ed is the energy of the isolated dot 
(which will also be referred to as the gate potential), c k 
(ct) is the destruction (creation) operator of an electron 
on the left (L) or right (R) lead with energy e k , and t k 
is the coupling amplitude between the dot and the lead 
level k. 

For a full specification of the model, one requires the 
spectral function which determines t k . In the applica- 
tions reported below, the leads will be described within 
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the wide band limit with a sharp cutoff at high and low 
energy values, such that the spectral density is given by 



Jl/r. (e) 
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L/R 
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fn the calculations reported below, unless otherwise 
noted, we use T L = T R = §, T = T L + T L , A = 5T,B 
. i — for. We also use a uniform discretization (Sex) to 
select the energies of the leads' states, and thus the cou- 
plings are given by: 



tk i e k) 
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Our aim is to use the ASM and MaxEnt method to cal- 
culate the measurable dc conductivity and its frequency 
dependence and compare the analytic continuation ap- 
proaches to exact results. For the model Hamiltonian 
given above, the current operator, /, is taken to repre- 
sent the current from the left lead to the dot (similarly, 
one can look at the right current), and is given by the 

change in occupancy of the left lead, / = e^rr-, where 

Nl = X) c 3c c fc; an( l a brief calculations shows that 
fcei 



d- S 



Cfc 
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With the above definition, the current-current corre- 
lation function can be calculated exactly in order to ob- 
tain the conductivity based on Kubo's linear response 
theoryP2l Let us assume that one knows the matrix which 
diagonalizes the Hamiltonian in Eq. (11) and denote it by 



U. By making a linear transformation for the operators 

Cfc = Y^UkaCa and c\ = J^^ak^a ( wnere fr° m now on 

a a 

the index k = will be refereed to that of the dot level) , 
we can rewrite the expression for C (t) in terms of the 
eigenfunctions of diagonal Hamiltonian: 



C(*) = irr{e-^e i **/e- < **/} 



r, 2 
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In the above, e n are the energy eigenvalues from the di- 
agonalized Hamiltonian H = ^ejcjcj, f (e n ) are the 

i=0 

Fermi-Dirac distribution at equilibrium, and 
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^t k t q (U m lU kn - U m lU n) (U m lU qn - U m lU 0n ) 
q,k 

(16) 

can be refereed to as the electron velocity matrix element 
in the the orthogonal basis, which obeys S„„ = for the 



same quantum state. The Fourier transform of Eq. (15) 



gives 
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where A nm — e n — e m and 5 n (x) is defined by 

5 n(*) = l 2 ,2 - (18) 
n x 2 + (hrf) 

Using Eq. ([5]), the experimental measured electrical 
conductivity can be calculated and is given by: 



S?KH}= lim °Y> 



' ■ run ^J mn S V ( A nm + fad) 

" Aran 



(19) 

where Af mn = / (e m ) — f (e n ). In practice, we use a 
finite value of ?y since the leads (baths) are described by 
a finite set of modes. We take hq to be larger than the 
spacing between adjacent bath modes and smaller than 
the band width. Its actual value is determined by starting 
from a large value of rj and gradually decreasing it until 
convergence occurs 

Formally, for an infinite number of bath modes one can 
take the limit r) -> 0+P 



^ sinh 
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+ hu), (20) 

and the dc conductivity is then obtain by taking uj to 0, 
which gives the well known Landauer formula for con- 
ductanc e 1 41 * 42 ! 
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IV. RESULTS 

The application of the ASM and MaxEnt method re- 
quires as input the imaginary time current-current corre- 
lation function, C{— ir) given by Eq. [8j In principle, the 
calculation of such correlation functions requires numer- 
ical techniques such as those based on path integration, 
where a variety of approaches can be extended to imagi- 
nary time P™*^ In the present case, due to the simplicity 
of the model Hamiltonian (cf., Eq. (JTT|)), instead of re- 
ferring to a specific numerical implementation of a path 
integration approach, we derive an exact expression for 
C(—ir). To preform the analytic continuation based on 
the ASM and MaxEnt method, we add artificial noise to 
the exact expression of C(—ir) to mimic the results of a 
Monte Carlo procedure, as described below. 

For a finite value of ij (see the discussion above), the 
imaginary time current-current correlation function is 
given by the exact expression: 
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where fl — j!m=lb* = iM a _ ^2 its \,r~ s pg n 

oo oo _ 

and .Ej (a;) = — J ^dt, E n (x) = — J ^-dt are the ex- 

l 

ponential integral function and the n th exponential in- 
tegral function, respectively. The above expression for 
C v (t) is well defined for < r < hf3 and reduces to 



C(- 



r;->-0+ 



nm 



(23) 



in the limit of ?y — > + . 

In Fig. [T] we plot the current-current correlation func- 
tion in imaginary time for two values of f3. The results 
were obtained for a finite number of bath modes (iVf,), 
where N b = 1000 for both the left (L) and right (R) 
leads. The results are shown for fvq — 0.04T. However, 
we find that within the accuracy of the numerical ana- 
lytic continuation, the inversion is not sensitive to the 
value of r\. In fact, on the scale of the plots shown in 
Fig. [I] different values of r\ (within a reasonable range) 
are indistinguishable. This implies that one can refer to 
Monte Carlo simulations to obtain the imaginary time 
data for a finite number of bath modes and ignore the 
role of rj in the simulations. This, of course, is not the 
case in real time, where for a finite system a K (u = 0) — > 
when j] = OP 




0.2 0.3 



0.4 



0.5 



Figure 1. Current-current correlation function in imaginary 
time (cf., Eq. (22} for /3 = 3/r (upper panel) and /3 = 10/T 
(lower panel). Results are shown for hrj — 0.04r and ed = 0. 
On the scale of the plots, it is difficult to differentiate between 
results obtained for different values of rj. 

Next, we applied both the ASM and MaxEnt method 
to invert the imaginary time data. We added Gaussian 
noise to the imaginary time current-current correlation 



function with a standard deviation of 1% at each data 
point. We discretized the imaginary time axis to N T — 
100 points and the frequency axis to N u — 1024 points. 
For the ASM approach, we averaged the solution over 
10 7 Monte Carlo sweeps. In Fig. ^ we plot the results 
for the frequency dependent electrical conductivity at a 
gate voltage = and 64 = T, and for two inverse 
temperatures of f3 = p and /3 
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Figure 2. Frequency dependent electrical conductivity o k {uj) 
obtained from the exact result from Eq. ( 19 \ (solid black 
curve), the ASM (dashed-dotted red curve) and the Max- 
Ent (dashed green curve) for different gate voltages (ed) and 
temperatures, ed is denoted in units of F, j3 in units of 1/r, 
and uj in units of V/h. Inset: Frequency dependent electrical 
conductivity a K {to) for the asymmetric coupling case Tl = \ 
and Y R = jj- y 

We find that the MaxEnt provides an overall better 
agreement with the exact results compared to the ASM. 
Both approaches provide reasonable description at low 
frequencies, and thus are accurate enough to determine 
the dc conductivity (see below), though MaxEnt is in 
general more accurate in this spectral regime. Both ap- 
proaches over-estimate the magnitude of the first peak 
observed at low frequencies. In general, MaxEnt does 
provide a better estimate of the position, height and 
width of the peak for the range of model parameters stud- 
ied in this work. The ASM also produces a sharp peak 
at higher frequencies which does not appear in the ex- 
act solution for a K (uj), while MaxEnt performs better in 
this regime. In fact, the high frequency behavior is quite 
tough for analytic continuation approaches. In D(uj) (not 
shown) we observe a sharp fall of the response as the fre- 
quency approaches the band cutoff B. This is translated 
to a small shoulder observed at hui = B in ct(uj). Such a 
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sharp change in the spectral behavior is challenging for 
the analytic continuation methods, and this explains the 
failure of the ASM. 

The inset of Fig. [2]shows a typical result for the asym- 
metric coupling case, where = \ and Tr = jk. This 
case is more difficult for analytic continuation since the 
overall magnitude of the conductivity is much smaller, 
due to the occurrence of destructive interference on the 
dot. MaxEnt provides reasonable results at all frequen- 
cies while the ASM fails markedly even at low frequencies. 
The fact that MaxEnt captures the behavior of the con- 
ductivity for the asymmetric case is encouraging. In fact, 
it suggests that the MaxEnt approach can be used for 
more elaborate situations: For example, when the cou- 
pling to the left lead is not proportional to that of the 
right lead and the simple Landauer type Meir-Wingreen 
formula is not adequate (see Eq. (9) in Ref. 




Figure 4. The dc conductivity a K (uu — 0) as a function of the 
gate voltage for temperature /3 — 3/T (upper panel) and 
P = 10/r (lower panel). 
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Figure 3. Same as Fig. [2] but for different values of the lead 
cutoff parameter, B, given in units of F. 



The high frequency behavior is also determined by the 
band structure of the leads. In Fig. [3] we plot the fre- 
quency dependent electrical conductivity for different val- 
ues of the band cutoff parameters, B defined in Eq. ( 12 1. 
Indeed, the MaxEnt method gives a well behaved spec- 
tral response while the ASM provides two peaks for all 
values of B. This is even more pronounced at small val- 
ues of B. where the width of the peaks in the ASM are 
very narrow, in contrast to the exact response which is 
smooth. Moreover, even the low frequency response of 
the ASM significantly deviates from the exact results at 
low band width, while within MaxEnt the accuracy at 
low frequency remains the same for all values of B, and 
the dc conductivity is only slightly affected by changing 
B. 



The dc conductivity as a function of the gate voltage 
and inverse temperature is shown in Fig. [4] and Fig. [5] 
respectively. The exact results given by Eq. (i"9| agree 
with the Landauer formula (not shown) HE For the cur- 
rent set of parameters, both approaches provide good re- 
sults for the conductivity at ui = 0. Again, the MaxEnt 
method seems to provide more accurate results near the 
resonance and is quantitative for most gate voltages and 
temperatures. 




Figure 5. The dc conductivity as a function of the inverse 
temperature for gate voltage Ed = T (upper panel) and e<j = 
(lower panel). 



V. CONCLUSIONS 

We have presented an analytic continuation approach 
based on Kubo's linear response theory to obtain the 
dc and ac components of the electrical conductivity in 
molecular junctions. Kubo's formulation requires the 
calculation of the current-current correlation function in 
real time, which is a difficult task for a many-body open 
quantum system due to the well-known dynamical sign 
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problem. The calculation of the corresponding imaginary 
time correlation, on the other hand, is a much simpler 
task, amenable to Monte Carlo techniques. 

Two approaches were adopted here to carry the ana- 
lytic continuation of the current-current correlation func- 
tion to real time: the ASM and the MaxEnt method. To 
assess the accuracy of these methods, we performed cal- 
culations for the resonant level model at a wide range 
of temperatures, gate voltages and frequencies. The nu- 
merical results were compared with an exact expression 
for the electrical conductivity, which is straightforward 
to obtain for this model. 

We find that MaxEnt is superior to the ASM for the 
entire range of frequencies and for different model pa- 
rameters. It provides an accurate description of the dc 
conductivity as well as a reasonable approximation for 



the ac component. Furthermore, MaxEnt captures in- 
terference effects on the dot resulting from breaking the 
symmetry in the couplings between the dot and the leads. 
The ASM fails in all these respects. 
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